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We investigate the cosmological dependence and the constraining power of large-scale galaxy cor¬ 
relations, including all redshift-distortions, wide-angle, lensing and gravitational potential effects on 
linear scales. We analyze the cosmological information present in the lensing convergence and in the 
gravitational potential terms describing the so-called “relativistic effects”, and we find that, while 
smaller than the information contained in intrinsic galaxy clustering, it is not negligible. We inves¬ 
tigate how neglecting them does bias cosmological measurements performed by future spectroscopic 
and photometric large-scale surveys such as SKA and Euclid. We perform a Fisher analysis using 
the CLASS code, modified to include scale-dependent galaxy bias and redshift-dependent magnifi¬ 
cation and evolution bias. Our results show that neglecting relativistic terms introduces an error in 
the forecasted precision in measuring cosmological parameters of the order of a few tens of percent, 
in particular when measuring the matter content of the Universe and primordial non-Gaussianity 
parameters. Therefore, we argue that radial correlations and integrated relativistic terms need to 
be taken into account when forecasting the constraining power of future large-scale number counts 
of galaxy surveys. 


I. INTRODUCTION 

At the beginning of the next decade, a variety of new galaxy surveys is expected to provide high precision data 
on galaxy clustering on very large scales. Planned and proposed experiments such as the Prime Focus Spectrograph 
(PFS) PP , DESI [21, SPHEREx [3], the Square Kilometre Array (SKA) [1], Euclid [5] and WFIRST [B] will probe 
cosmological volumes an order of magnitude larger than what has been provided by current surveys, so there is the 
need of a precise theoretical modeling of galaxy clustering. 

Analyses of power spectrum and correlation function have been performed mainly using the so-called Kaiser for¬ 
malism 00, that assumes (other than linearity) the plane-parallel approximation, and it also sets the two galaxies of 
the pair at the same redshift. Classical analyses are also performed using Newtonian approximations, and the validity 
of this approach breaks down when probing scales approaching the horizon, where quantities have to be defined in 
a gauge-invariant way consistently including perturbations of observed redshift and volume. Theoretical models for 
galaxy clustering on large scales were developed and tested including wide angle corrections to the Kaiser approxi¬ 
mation and radial effects [SHm, and, more recently, including all effects due to the fact that we observe galaxies on 
the physical, perturbed past light cone and not on some fictitious Friedman background [20H32] , 

Future galaxy surveys are an excellent opportunity to advance our understanding of cosmology in general and Large 
Scale Structure (LSS) in particular, e.g. via constraining the Dark Energy equation of state or by testing General 
Relativity jT£T, 31] , Furthermore, joint analysis with CMB experiments |351 30 will bring better constraints on the 
primordial Universe. Inflation models can be constrained studying, e.g., the primordial power spectrum and features 
induced by primordial non-Gaussianity. The latter affect the galaxy correlation function on very large scales. The 
interesting results of EUlH] show that primordial non-Gaussianity are detectable also via the two-point correlation 
function (or its Fourier transform, the power spectrum) since they enter as a scale-dependence of the large-scale linear 
galaxy bias. 

The main goal of the present paper is to study the cosmological dependence of relativistic effects and the error 
in parameter estimation introduced when those terms are neglected. We will also investigate which term is mainly 
responsible for the results obtained in [39 ] 140 ] . and we extend their previous analysis to other cosmological parameters. 
We consider the SKA and Euclid surveys, consistently including galaxy bias, magnification bias and evolution bias. 
We study the number counts angular power spectra at first order in perturbation theory and perform a Fisher matrix 
analysis to estimate cosmological parameters. Our aim is not to forecast the most stringent and precise error bars for 
cosmological parameters, but to compare the standard Newtonian analysis to the fully relativistic one. For this we 
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assume rather conservative specifications. This yields relatively large error ellipses, but assures that our analysis is 
not affected by, e.g., the treatment of non-linear scales or by too optimistic survey specifications. A detailed analysis 
of the goodness of the constraints per se is left as a future work including forecasts based on Markov chain Monte 
Carlo methods. 

The paper is organized as follows. In Section [II] we introduce the formalism used throughout the paper and 
in Section m we present the surveys specifications used to model number density as a function of redshift, the 
galaxy bias, the magnification, the evolution bias and discuss our Fisher analysis. In Section |IV| we investigate how 
relativistic effects can be used to constrain cosmological models and what is their effect on parameter estimation, and 
in Section [V] we show that relativistic effects contain cosmological information, and so they can in principle be used to 
constrain cosmological models. Our results on how neglecting relativistic effects can modify the predicted cosmological 
measurements of future surveys are summarized in Section |VH and we discuss our conclusions in Section [VlI| Appendix 
[A| describes our implementation of primordial non-Gaussianity into the Class code. 


II. GALAXY CORRELATIONS IN GENERAL RELATIVITY 

We cannot observe galaxies in real physical space at constant time. We observe them on our background light cone. 
Actually, for each galaxy we observe its direction n in the sky as well as its redshift z. Taking into account that 
this redshift is perturbed by peculiar velocities is the well known “redshift space distortion” originally pointed out 
by Kaiser [7j. However, also the gravitational potential induces an apparent redshift modification and the direction 
into which we see the galaxy may not be the one into which the photons have been emitted. Taking into account 
all these relativistic effects to first order, one can define the observed galaxy over-density A Q b s (n, 2 ) at fixed observed 
redshift and into a given observed direction n. This is an observable and therefore gauge-invariant quantity which 
involves not only the density fluctuation and the velocity but also the Bardeen potentials <f> and 'F. The fact that this 
expression is more complicated than the simple Kaiser formalism is not only a difficulty, but as we shall show, also 
an opportunity. It means that the number counts contain information not only on the density but also, e.g., about 
lensing, which is relevant and does affect cosmological parameter estimation. 

Real and Redshift Space in GR 



Hubble expansion + 



FIG. 1: Left Panel: We observe galaxy directions and redshifts on our perturbed past light cone. Right Panel: Effect on the 
apparent position of a galaxy: peculiar velocity, gravitational potential and lensing modify the real position and velocity (in 
the center of the image, circled), to the observed one . 

The total overdensity A 0 b s (n, z) is not anymore just the local density field, but it is comprised also by integrated 
terms that correspond to lensing, time-delay and integrated Sachs-Wolf effects, due to the impact of cosmological 
perturbations on the path of photons. We can then express the total observed over-density as: 


A ob s(n, z) = A 5 (n, z) + A rsd (n, z) + A v (n, z) + A K (n, z) + A pot (n, z ), (1) 

where 5 refers to the overdensity in the comoving gauge, rsd and v are peculiar veolcity (redshift space distortions in 
the Kaiser approximation) and Doppler effects, k contains lensing convergence and pot incorporates local and non local 
terms depending on Bardeen potentials and their temporal derivatives. The definition of these in terms of standard 
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linear perturbation variables are given in Appendix [Aj There we also present the transfer functions of these terms, 
including galaxy bias b(z , k), magnification bias s(z ) (in the literature also the quantity Q{z) = 5s(z)/2 is sometimes 
used) and evolution bias f eV o{z) (sometimes called b e (z)). 

In Figure [T] we illustrate the apparent modifications to the position of a galaxy due to all the effects mentioned; 
in the “standard” case the only effect present is the one due to the Kaiser approximation for peculiar velocities that 
causes a spherical distribution in comoving coordinates space, to be observed as squashed along the line of sight 
(see e.g. Figure 1 of [8]). In the relativistic treatment there are also contributions from lensing convergence k and 
gravitational potential terms, that cause an apparent modification to the angular and radial position, respectively. 
The convergence n is normally anti-correlated with the overdensity 6 but the magnification bias which enters with the 
opposite sign can dominate it, especially at high redshift, see, e.g. Figure 12 of [41]. Because of this, a distribution 
of galaxies may appear more or less compact. Finally, the gravitational potential terms directly involve Bardeen 
potentials, their time derivatives and their integrals along the line of sight (time-delay and integrated Sachs-Wolf 
effects). 

These correlations have been investigated recently in [HI [231 EH [29], and it has been shown that, for radial 
correlations with a large A z, relativistic terms can be a significant contribution to the correlation m- 

Since A 0 b s (n,2) is a function of direction and redshift, its power spectra are defined for two redshifts 2 and z ', 
Cf{z. z'), see [ 21 ] ■ In Appendix |a| we define the transfer function A^(k) of the number counts given in Equation |lj. 
This computes the contribution to A 0 b s from wave number k estimated in a bin with mean redshift Zi and window 
function R). From the transfer function, the angular power spectra are determined by integration over the primordial 
curvature spectrum [29] : 


cy =4nf ^ Af* (fc)A^ (k)Vn(k) , (2) 

where Viz(k) is the primordial power spectrum of curvature perturbations. The superscript Wj indicates that these 
transfer already include an integral over the window function IT). More details are given in Appendix |A| 

We note that in the standard analysis of power spectra P(k) it is customary to use only auto-correlations within a 
certain redshift bin of a galaxy survey, and not the cross-bin correlations as redshift bins are treated as independent. 
In the following we show that it is important to be able to model consistently cross-bin correlations, since they are 
dominated by the non-negligible lensing convergence n. This is very interesting since k is the (spherical) Laplacian 
of the lensing potential ip, hence their spectra are related by Cjp = [£(£ + 1 )] 2 Cf. Therefore, if galaxy surveys are 
sensitive to n, we can in principle use them to determine the lensing potential, usually inferred by shear measurements. 


III. FUTURE GALAXY SURVEYS 

To investigate the observational consequences of including relativistic corrections we consider two examples of future 
surveys: the HI galaxy survey planned with the SKA m as an example of a spectroscopic survey, and the Euclid 
photometric survey [5] [33]. We use specifications for SKA which are consistent with the 5/iJy cut of [13, and the 
photometric specifications for the Euclid survey. In Figure [2] [3] we plot the predicted number of galaxies, galaxy bias 
and magnification bias that we will use for our forecasts. 


A. Survey specifications 


1. Galaxy distribution 


The number of galaxies per steradian and per unit redshift for SKA [48] is given by: 


dN 

dzdil 



2 

10 Cl ^ C2 e" C3Z 


(3) 


we use Ci = 6.7767, c 2 = 2.1757, C3 = 6.6874; the survey covers 30,000 deg 2 , i.e. / s k y = 0.73. We divide the survey 
into 5 bins of equal depth Az = 0.39 in the range 0.05 < z <2. 

For the Euclid photometric survey we use: 


dN 

dzdfl 


= 3.5 x 10 8 ,z 2 exp 



0 < z < 2, 


(4) 
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FIG. 2: Redshift distribution, normalized to unity, (left panel) and bias (right panel) for the SKA and Euclid surveys considered 
in this paper. 


where Zo = z mea n/l-412 and the median redshift is z mean = 0.9; f s k y = 0.375. This number density is consistent with 
the official specification for the density of galaxies, d = 30arcmin -2 . 

In both cases, these specifications will give a total number of sources detected of ~ 10 9 . 


2. Galaxy bias 

The galaxy bias for SKA is approximated as j48j : 

b G (z) = c 4 e C5Z , (5) 


with a = 0.5887, C 5 = 0.8137. 

As for Euclid, we model the galaxy bias as: 


b(z) = Vl + z . 


( 6 ) 


These fits are conveniently taken into account in the implementation of primordial non-Gaussianities discussed in 
Appendix 


3. Evolution bias 

The rsd and pot terms also include evolution bias, f evo (z), (see e.g. [251 [Ml I5U] ). This is due to the fact that new 
galaxies form so that the true number density of galaxies does not simply scale like a -3 where a is the cosmological 
scale factor. The evolution bias is given by: 


/e vo (z) 


— — In (a 3 dN ( z ’ L > Llim 2 
aH dr \ dzdfl 


(7) 


where r is the conformal time and H is the Hubble parameter. Here dN(z , L > An m ) /dz / d£l indicates the true number 
density of galaxies (not necessarily observed) per redshift and per solid angle present in the Universe above a certain 
threshold L\\ m in luminosity. It can be estimated from the luminosity function |22| . However, given the uncertainties 
in the modelling of galaxy evolution , for simplicity as in [39] we assume that the observed dN(z)/dz/dfl in Figure [ 2 ] 
still gives a good approximation to estimate f eV o( z )- This comes from the fact that the evolution bias only appears 
in subleading terms (wide-angle velocity and pot) so that uncertainties in its modeling do not significantly affect our 
results. 
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4- Magnification bias 

For SKA we interpolate the results of [39] to 5/iJy: 

s(z) = c 6 + C 7 e~ CsZ (8) 

with c 6 = 0.9329, C 7 = —1.5621, c 8 = 2.4377, and use a constant s(z{) within each bin centered at z ir 
We use the magnification bias for the photometric Euclid survey computed in [84] . given by the fit: 

s(z) = s 0 + Siz + s 2 z 2 + s 3 z 3 (9) 

with So = 0.1194, Si = 0.2122, S 2 = —0.0671 and s 3 = 0.1031. This result is consistent also with other analyses, 
e.g. US]- In Figure [3] we plot the magnification bias used in this work. We modified the CLASS code to introduce a 
magnification bias s[Zi) depending on the central redshift Zi of each bin. 




FIG. 3: Magnification bias for the SKA (left panel) and Euclid (right panel) surveys considered in this paper. 

Cosmic magnification changes the relative number of sources detected at a given redshift and fixed magnitude limit, 
and in this way lensing can modify the effect of cosmological parameter variations on clustering (see e.g. CM): 

n 0 b s {z) = n g {z)[ 1 + S + (5s(z) - 2 )k] , (10) 

where n 0 \, s ,n g are the observed and intrinsic number of sources, respectively, s is the magnification bias and k is the 
convergence. We stress that this effect modifies the effect of primordial non-Gaussianity also, due to changes in the 
relative redshift distribution of observed sources, and so lensing effects are important for measurements of /nlj a fact 
that is not often taken into consideration. 


B. Fisher matrix analysis 

In order to predict the constraining power of a survey one often uses a Fisher matrix analysis [±I il • The Fisher 
matrix, F a p, indicates how well a given set of cosmological parameters ($ a ) under consideration can be measured. 
The bigger the Fisher matrix elements, the better the parameters can be determined by the given survey. If the errors 
(in the parameters) would be Gaussian, the ellipsoid in parameter space given by: 


(«=i 


X! F c*p(tia ~ $aWp ~ & p) < A ra y 2 


aB= 1 


where the Inverse Regularized Gamma function gives the difference to the minimum y 2 as Aiy 2 = 1 and A 2 X 2 = 2.30 
for n = 1,2 degrees of freedom, respectively. This determines the 1 <t error bars or 68 % confidence region of parameter 
space around some fiducial values (i9 a ). In reality errors are not Gaussian but the Fisher ellipses usually still give a 
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reasonable indication of the precision with which a given parameter can be determined. In our situation, where we 
measure the C\° with a certain precision ac e , the Fisher matrix is given by: 


p _ Y= dCJj^dCf^o 

“ 0 ^ ^ dd a dda Cl{{ij) ' {pq)V 


( 11 ) 


where the derivative is evaluated at fiducial values and ac e are errors in the power spectra. For Gaussian fluctua¬ 
tions they can be determined by Wick’s theorem and are given by (see e.g. [29]): 


^y£ ^y £ i {y £ {y t 

(2£ + l)/ sky 


iq)/=*(jp) 


a c, 


<[(U).(pq)] 


( 12 ) 


where Cp denote the observed correlation multipoles which include shot noise errors: 

6 ' = c ' + JW^n- 


(13) 


and dN(zi)/dfl is the average number of sources per steradian within the bin z,. Even with primordial non Gaussianity 
this is a good approximation for the error in the power spectra. Note that we sum over the matrix indices ( ij ) with 
i < j and ( pq ) with p < q which run from 1 to the number of bins. 

To perform the Fisher analysis, we parameterize our cosmology using the following parameters: 


P — {W 0 , W a , 77. s , G( s , flcdm; Hbj /nL; ^Ng} • 


(14) 


We have set spatial curvature K = 0; h parameterizes the present Hubble parameter, H$ = M00km/s/Mpc; {wo,w a } 
parameterize the dark energy equation of state, w = wq + w a { 1 — a); {n s ,a s } are the primordial power spectral 
index and its running (see definition in the next section); {f2 cc j m , n b } are the density parameters of cold dark matter 
and of baryons respectively; {/nl,?ing} are the primordial non-Gaussianity parameter and its spectral tilt (see 
definition in the next section). We choose the fiducial values u>o = —1, w a = 0.05, d s = 0.01, ting =0.1, /nl = 1, 
h = 0.67, H cdm = 0.25 and S2 b = 0.05. 

In doing a real data analysis, one should use the range 2 < i < fcmaxX(^), but this requires involved algorithms to 
carefully exclude non-linear scales [29] and to compute the beginning of the quasi-linear regime when cross-correlating 
different redshift bins. Here we are mostly interested in a first assessment of the importance of relativistic effects as 
compared to the standard analysis, so to simplify our study, we assume a conservative ^ max =200 at all redshifts. 
Even if some of the relativistic effects are more sensitive to low multipoles, such a conservative cut on non-linear scale 
does not bias our constraints towards an ad-hoc enhancement of relativistic effects compared the Newtonian case. In 
fact some of the most relevant terms (e.g. lensing) also affect high multipoles, see Figures |4j [5] 


IV. ESTIMATING COSMOLOGICAL PARAMETERS WITH RELATIVISTIC GALAXY 

CORRELATIONS 

The main goal of this paper is to investigate the impact of relativistic corrections on galaxy clustering analyses and 
their constraining power by looking at the predicted Ce(z,z r ) from future galaxy surveys. We focus on constraints 
on dark energy and on the parameters defining the initial perturbations. Note that we are studying the impact of 
relativistic terms relative to the Newtonian ones and we do not perform a forecast of survey capabilities; a detailed 
analysis of survey specifications and observational strategies together with an MCMC analysis including the nuisance 
parameters of the survey is required to correctly assess the proper forecasted measurements for specific instruments. 
This is not the focus of this work. 

In this paper we compare constraints using a traditional analysis (but including wide-angle and local relativistic 
corrections; see m for more details) with a fully relativistic one. We want to study how much information is present 
in the relativistic terms which have been neglected in the past. In order to do this, we compare two Fisher matrix 
analyses: the Newtonian case, referred to as “nwt” where we include intrinsic clustering plus peculiar velocity effects 
(in the formalism of Equation [lj this corresponds to the first two terms), using only auto-correlation functions within 
the redshift bins used, and the fully relativistic case, referred to as “rel” where we include all the terms and also 
cross-bin correlations. 

In Figure [4j [5] we compare the angular power spectra in the Newtonian and relativistic case for the SKA (spec¬ 
troscopic) and Euclid (photometric), respectively. In the plots we refer to < zaZb > to indicate the correlations 
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Ce(zA, zb)- In the case of spectroscopic surveys we assume top hat bins, while we use Gaussian bins in the case 
of photometric surveys; more details on the modeling and computation of Ce(zA, Zb) are given in Appendix |A} In 
both cases, small differences are visible in the auto-correlations, with larger differences at small i\ on the cross-bin 
correlations the differences are very large, because these correlations are dominated by the lensing convergence term. 
In particular, non adjacent bin cross-correlations are nearly zero in the newtonian case for spectroscopic surveys, 
as expected. We note that oscillations in e.g. the {z\Z±) case plotted here are not physical but due to numerical 
precision for such low (negligible) amplitudes. This justifies the standard assumption that in spectroscopic surveys 
the cross-bin correlations contain virtually no cosmological information, but, as firstly noted in m , we stress the fact 
that in the relativistic case, gravitational potential and lensing convergence effects introduce radial cross-bin large 
correlations. Nearest cross-bin correlations, e.g. ( 2425 ) from Euclid are in general larger compared to the SKA ones in 
the Newtonian case, because photometric redshift uncertainties assumed for Euclid cause an effective cross-bin overlap 
(see e.g. jHEj)- This means that local Newtonian term are still non-negligible compared to the integrated relativistic 
ones. We will investigate the cosmological information present in these correlations in the rest of this Section. 





FIG. 4: Power spectra in the Newtonian and Relativistic cases for different Ci(za, zb) for the spectroscopic SKA galaxy 
survey. Differences in the auto-correlations are small and at very low £, while cross-bin correlations are dominated by integrated 
relativistic terms. 


A. Dynamical Dark Energy 

The first step towards understanding the nature of dark energy is to clarify whether it is a simple cosmological 
constant or whether it originates from a dynamical degree of freedom, i.e. it varies with time. Dynamical dark energy 
can be distinguished from the cosmological constant by considering the evolution of the equation of state of dark 
energy, w = p , where p and g are the pressure and energy density of the “dark energy fluid”, respectively. For a 
cosmological constant model, w = — 1 , while for dynamical models w = w(a), where a is normalized to unity today, 
a 0 = 1 . 
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FIG. 5: Power spectra in the Newtonian and Relativistic cases for different Ce(zA, zb) for the photometric Euclid survey. 
Differences in the auto-correlations are small and at very low l, while cross-bin correlations are dominated by integrated 
relativistic terms. 


To evaluate the potential to constrain the dynamics of different models of dark energy with galaxy surveys, we 
adopt the following simple phenomenological parameterization for the dark energy equation-of-state (EoS) w [891190] : 

w(a) = Wq + w a (l — a). (15) 

In Figure [6] we show the constraints on these dynamical dark energy parameters for the SKA and Euclid, marginal¬ 
ized over other cosmological parameters (here and throughout the entire paper we will show l-er constraints). It can 
be seen that, as expected, differences in the estimation of the errors are larger in the photometric case. We will 
quantify the error on the forecasted errors in the determination of dynamical dark energy parameters in Section |VI| 
For spectroscopic surveys, including lensing, and especially when magnification bias s(z) is consistently included, 
actually degrades the precision with which (wo,w a ) can be determined by this observation. 


B. Initial conditions from Inflation 

We study how large-scale correlations can be used to test the initial conditions of scalar perturbations. In particular, 
we focus on the spectral index, n s of the power spectrum and its running, a s as well as the primordial non-Gaussianity 
parameter /nl and its spectral tilt n^G- 


1. The primordial Power spectrum 

Inflation predicts a spectrum of primordial curvature and density perturbations that are the seeds of cosmic structure 
which forms as the universe expands and cools. This initial spectrum, and in particular its spectral index n s is a key 
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FIG. 6: Constraints on dynamical dark energy parameters {wo,w a } from SKA (spectroscopic, left) and Euclid (photometric, 
right). 


observable for understanding inflation. The spectrum is defined by: 

k 3 (TZ( k)ft(k')) = d(k - k ')P K (k ). 


(16) 


The Dirac-delta is a consequence of statistical homogeneity and because of statistical isotropy Pr only depends on 
the modulus of k. The spectral index n s is given by: 


n a - 1 = 


dlnP-R 


dlnfc 


fc=fc» 


(17) 


where /c* denotes the pivot scale which is usually chosen to be a scale to which a given experiment is most sensitive. 
Given this, and to simplify the comparison with CMB experiments [92], we choose fc* = 0.05/Mpc. The quantity n s — 1 
is also called the spectral tilt. Inflation predicts a nearly scale-invariant spectrum of primordial density fluctuations, 
i.e. n s ss 1 and a nearly scale independent tilt, i.e. n s nearly independent of the pivot scale. We characterize a 
possible scale dependence of n s by: 


dn s 


cg = 


din (fc) 


(18) 


k—k^ 


the quantity a s is called the “running” of the spectral index. 

Figure [7] shows the constraints on power spectrum spectral index parameters, marginalized over the other cosmo¬ 
logical parameters. The effects of integrated terms is in this case smaller, and in the case of spectroscopic surveys, 
neglecting relativistic effects will induce only a ~ 5% error in the estimated precision, as we will see in Section VI It is 
however worth noting that the results here are obtained by marginalizing over the remaining cosmological parameters, 


while in Section VI we fix them. In the photometric case, on the other hand, relativistic effects contribute significantly 
to the information content of galaxy clustering over the entire survey. 


2. Primordial non-Gaussianity 

An important goal for forthcoming cosmological experiments is to test whether initial conditions of cosmological 
perturbations deviate from Gaussianity. This is very interesting as standard slow roll single field inflation predicts 
very small non-Gaussanities [91]. This can be tested using CMB data ([36, 55, 56] and references therein) or using 
galaxy catalogs [32l|37H40l[46l[57H67]. Constraining primordial non-Gaussianity offers a powerful test of the generation 
mechanism of cosmological perturbations in the early universe. While standard single-field models of slow-roll inflation 
lead to small departures from Gaussianity, non-standard scenarios (such as e.g. multi-field inflation) allow for a larger 
level of non-Gaussianity (see e.g. [Ml El]). A widely used parameterization of primordial non-Gaussianity is the /nl 
parameterization, that includes a quadratic correction to the potential: 

$NG = <t> + /nl (< t >‘ 1 - (/> 2 )) , 


( 19 ) 
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FIG. 7: Comparison of the error ellipses for the spectral index and its running in the Newtonian (dashed) and in the relativistic 
analyses (solid) for SKA (spectroscopic, left panel) and Euclid (photometric, right panel). 




where <&ng denotes the Bardeen potential, which, on sub-Hubble scales reduces to the usual Newtonian gravitational 
potential. Here <j) is a Gaussian random field, and the second term, when /nl is not zero, gives the deviation from 
Gaussianity. In this paper we refer to the so-called “local type” /nl- We use the LSS convention (as opposed to the 
CMB one, where ~ 1-3 /nl 10 see [BO] and Appendix |a| . 

One method for constraining non-Gaussianity from large-scale structure surveys exploits the fact that a positive /nl 
corresponds to positive skewness of the density probability distribution, and hence an increased number of massive 
objects [371 EH EZI ■ In particular, /nl introduces a scale-dependent modification of the large-scale halo bias. The 
difference from the usual Gaussian bias is given by [33]: 


= NM - i]/NL(t)4 v ^°"g (z) , (20) 

where bc(z) is the usual bias calculated assuming Gaussian initial conditions, which we assume to be scale-independent, 
D(z) is the linear growth factor, flo m is the matter density parameter today and 5 ec is the critical value of the matter 
overdensity for ellipsoidal collapse. We choose to approximate it by the spherical collapse value S sc ~ 1.68 lB31 . 

While generally the parameter /nl is assumed to be scale-independent, some inflationary models predicts a scale- 
dependent /nl> see e.g. IMIS] for more details. The running parameter ting is expected to be < 0(1) in most 
inflationary models. Its signatures in the CMB and LSS have been discussed in the literature, see e.g. [73H3T] . We 
consider the scenario that isolates the multi-field effects as in J33, S3], for which the effective /nl(^’) becomes: 

/ u \ 

/nl(*0 = /nl ( T -) > (21) 

\ fc *,NG / 

where we choose the pivot scale /c^ng = 0.04/Mpc to easily compare analyses of the CMB [33]. Effects of primordial 
non-Gaussianity on the bias are relevant on large scales due to the factor k~ 2 like the relativistic effects [SB], so we 
may encounter a certain degree of degeneracy. However, as their angular dependence is different, they can in principle 
be distinguished in a multipole expansion [26] . 

Figure [ 8 ] shows the constraints on primordial non-Gaussianity parameters, marginalized over other cosmological 
parameters. We can see that in this case neglecting large-scale effects has a large impact, confirming, with a different 
analysis, the findings of mm- We expect that in general lensing convergence is the dominant term in the radial 
cross-bin correlations (see 123121]), but since lensing is a tracer of the matter distribution and insensitive to galaxy 
bias, it is usually not taken into account when estimating non-Gaussianity from LSS. However, cosmic magnification 
effects, as from see Equation [lOj modify the effect of /nl on clustering. 

The variation in the ellipses is again larger in the photometric case, but it is still of ~ 20% for spectroscopic surveys. 
We stress again that our results here are not a prediction of the final constraining power of future surveys. Several 
differences between our analysis and the ones of e.g. [33] [63] give us different predicted constraints for the SKA (e.g. 
the number of bins, number of parameters over which we marginalize, inclusion of running of /nl)- Given the precision 
in the measurements of non-Gaussianity parameters that is forecasted for future galaxy surveys when performing an 
optimized analysis (3, 63) [65], a ~ 20% effect on the estimate of the error can be important for drawing theoretical 
conclusions from observations (see e.g. [BBj). 
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It is also worth noting that, as can be seen from Figure [8j including relativistic effects changes the degeneracy 
between /nl and ting- 




FIG. 8: The same as Figs. [6] and [7] for primordial non-Gaussianity parameters. 


V. COSMOLOGICAL INFORMATION IN RELATIVISTIC CORRECTIONS 

In this Section we investigate the cosmological information present in the relativistic corrections, given by the 
lensing and gravitational terms of Equation ([T]) . 

As stated in Section [TT| lensing and gravitational potential contributions modify the observed correlations, acting 
as corrections to the intrinsic clustering. Here we study the information contained in these relativistic corrections, 
by isolating them and performing a Fisher analysis using signal coming only from these effects (i.e. excluding galaxy 
clustering signal). 

In Figure [9] we show the constraints on cosmological parameters obtained by isolating the relativistic contributions. 
We marginalize over the cosmological parameters not considered in the given ellipse. Clearly, these terms depend on 
cosmological parameters, so they can in principle be used for testing cosmological models, if they can be isolated; but 
more importantly, neglecting them in an analysis causes a bias of the results. 

The constraining power of the integrated terms is of course smaller than the analysis that include galaxy clustering, 
but we can see how there clearly is some cosmological information in these correlations, and this explains the change 
it the ellipses in Section |TV| 

It is interesting to note the nearly perfect degeneracy in the {u>o,u; a } plane. Lensing is very sensitive to the 
combination wq + which it constrains to be —1 ± 0.03. 

These results also suggest that futuristic surveys or improved analyses (e.g. using the so called multi-tracer tech¬ 
nique, or specifically planned observations targeting galaxy samples with a very large magnification bias) could provide 
a complementary test for cosmological models. In particular, integrated terms may contain additional information 
on the model of gravity; an analysis of this for the f(R) model can be found in [25], and we will present a general 
investigation in a follow-up paper. 


VI. SUMMARY OF RESULTS 

To quantify and summarize our results, we construct Figures of Merit (FoM) describing parameters of cosmological 
models. These are given by the inverse of the area of the error ellipses shown in the previous figures, see e.g. [2S] for 
details. To properly quantify the error in the forecasted precision in measuring cosmological parameters introduced 
when neglecting relativistic effects, we then define the ratio: 

ta_ tv rrel ^ 

rF ° M = ’ F ° M = M ^}) -1 ])' 5 > ( 22 ) 

where { A , B} = {/nl; HigJi { w 0 ; w a }, {n s , a s }, {f2 c d m , ^b} describe the model parameters we want to test, and Fs^ ,b} 
is the sub-matrix of the Fisher matrix of the parameters {A, B} . This corresponds to fixing other parameters; clearly, 
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FIG. 9: Constraints on {n a , a 3 }, {wo, w a } from the integrated terms in Equation 0 for SKA (black) and Euclid (blue). 


the larger the FoM the better the given parameters are determined by the experiment. We compute the FoM for 
the spectroscopic and photometric example surveys (SKA and Euclid, respectively). In Figure [To] we show the errors 
in the in the predicted constraining power when performing a full relativistic analysis compared to the standard 
Newtonian analysis for all parameter subsets. 

1.5 


1.4 


1.3 

2 

o 


1.1 


1.0 

FIG. 10: Modifications in the Figure of Merit for different sets of cosmological parameters when going from the Newtonian to 
the relativistic analysis, for the photometric Euclid and the spectroscopic SKA galaxy surveys. 

We can see that in most cases there is a significant difference in the FoM for cosmological models when using a full 
relativistic analysis rather than a standard Newtonian one. This shows that a precise calculation of the FoM for future 
surveys requires a proper relativistic formulation of galaxy clustering. The largest effects are on the estimate of the 
matter content of the Universe and parameters describing primordial non-Gaussanity. As described in the previous 
Sections, effects are larger for photometric surveys, due to the fact that photometric bins will cause an increased 


SKA (spectroscopic) 


° {fNL> n NGl 

• {w 0 ,w a } 
x K,a s } 

D {^CDM^bl 


Euclid (photometric) 
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importance of the integrated radial terms coming from lensing and gravitational potential with respect to the intrinsic 
galaxy clustering signal. 


VII. CONCLUSIONS 


In this paper we investigate the effect of including relativistic corrections to galaxy clustering on the constraining 
power of future galaxy surveys such as SKA and Euclid, by using the directly observable angular power spectra. We 
study how cosmological parameter estimation is affected by relativistic terms by comparing the fully relativistic result 
with a standard Newtonian analysis. 

We also show that cosmic magnification together with other relativistic terms, via their effective “redistribution” of 
galaxies in redshift space, introduce a bias on the constraints on cosmological, including non-Gaussianity parameters, 
and so it will need to be taken into account when forecasting measurements of e.g. /nl- 

We analyze separately the integrated relativistic terms to determine their individual contributions to the signal and 
we show that they contain cosmological information. While the constraining power of the relativistic terms is smaller 
than the one of the standard density and redshift space distortion terms, neglecting them leads to a considerable bias 
in the estimated precision in the measurements of cosmological parameters describing the models for dark energy 
and for the initial perturbations. We quantify our findings in Section VI and we show that the error in neglecting 
relativistic terms is on the order of tens of percent, and of ~ 20 — 40% in most cases considered in this work. Future 
surveys aiming to measure cosmological parameters with high precision and to set strict constraints on theoretical 
models will then need to take properly into account radial correlations and relativistic terms in order to not bias their 
results. 

We don’t perform a detailed analysis of the constraining power of the surveys, as that is beyond the scope of this 
work; we expect the FoM to improve with the number of bins, that we keep to 5, in order to speed up computations 
and to have an easier understanding of the different contributions to the final results. Even though in slimmer bins the 
shot noise is increased due to the smaller number of galaxies, also the signal is increased since there is less “smearing 
out” of the dominant fluctuations on smaller scales so that the noise increase can be compensated (see [55] for details). 
It is an interesting future project to see how much can be gained by significantly increasing the number of bins. 

Finally we want to note that the angular power spectrum approach of projecting sources on a single value of z 
within the redshift bin misses the radial information in it. 
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Appendix A: Implementation of primordial non-Gaussianity effects in Class 

In this Appendix we summarize our modifications to the Class code ram to include effects of primordial non- 
Gaussianity effects, and show the number counts transfer functions computed in the code. Before we enter into this 
we want to clarify the definitions given at the beginning of the paper. 
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1. Definitions and notation 


The relativistic expressions for the number counts in linear perturbation theory is given in Equation |l]), 
A obs (n : z) = A s (n, z) + A rsd (n, z) + A v (n, z) + A K (n, z) + A pot (n, z ), 

where [ 2 TJ [ 22 j 

A^(n, z) = b(z,k) 6 co (r(z)n,T(z)) 


A rsd (n, z) = —-d r (V-n) 
U{z) 


A v (n, z) = 


W 2 — 5s(z) 


H 2 


rH 


+ 5 s(z) - fevo(z) 


(V • n) + [3H -/evo(*)] A-^V ■ V) 


A K (n, z) = (2-5 s(z))k = - 


(2-5 a(z)) r r(z) , r(z) — r 


dr ' ' A 2 (4> + W) 
ryz)r 


-y 0—5s 

A pot (n, 2 ) = (5s(z) - 2)<h + 'T + 'H~ 1 & + + ——- 1 - 5s(z) — f evo (z) 


H 2 rTL 


rr(z) 

T + / dr{& + T') 
Jo 


2-5 s f r{z) , T •. 

+ ^wi dr{ * + ^- 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A 6 ) 


Here V is the peculiar velocity, $ and T are the Bardeen potentials, 5 co is the density contrast in comoving gauge and 
% = aH is the conformal Hubble parameter. All quantities which are not integrated are evaluated at conformal time 
t(z) and at position r(z) n = (tq — r(z))n. Here r(z) is the conformal distance on the light cone, r(z) = t 0 — t(z). A 
prime indicates a derivative w.r.t. conformal time. The term A _1 (V • V) is the velocity potential. It comes from the 
non-local requirement that the overdensity (5 co which is multiplied with the bias factor b(z ) be the comoving one. In 
A:—space this becomes simply —V(k)/k. 

Expanding these expressions in spherical harmonics and defining by Sx(k,r) the solution of linearized Einstein 
equations for the variable X with initial condition $(/c,Ti n ) = 4/(fc, Tj n ) = 1 and >I> , (fc,Ti n ) = ^'(k,Ti n ) = 0, where 
Ti n is such that the corresponding mode is super horizon, i.e. kr ln C 1, we find after integration over a normalized 
window function IT) (z) centered at Zi that the transfer function A l ( Vi given in Equation ([ 2 ]) is the sum of the following 
terms (see [301 f° r more details): 


a Deni 


a Lerii 

aY 1 - 


aY 2 - 


aY 3 - = 


G 2 i 

I 

G3i 


/ drWi [b G (z ) + A b(z, k)} S s jn 
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14 —7777 T 7 - 7—77 T 5s — f e 
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r T ° Sa 

J dr Wi (/evo - 3 aH) -A- j e 

Sef; 


Se ., 

-Y J * 


[ dr Wi Sq, ]i 

/ 0 


(■TO 


dr Wi 


H’ 
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Jo 
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(A7) 
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Here we have also allowed for a scale dependent bias A b(z,k) which can be added to the Gaussian bias bc(z). We 
shall introduce this in the next subsection. Furthermore, also the functions s(z) and f ev0 (z) depend on redshift as 
described in the main text. 

As in Ref. [3D], the different contributions correspond to the density in comoving gauge 5 co (“Den”), lensing 
convergence k (“Len”), Doppler (“V1”-“V2”), redshift-space distortions in the Kaiser approximation (“V3”) and 
terms depending on gravitational potential (“G1”-“G5”), respectively. We use conformal time r, the proper Hubble 
parameter H and the scale factor a. We have omitted the arguments k for the transfer functions A, (r, k) for the 
source functions Sx, rk for the Bessel functions jt, and r for selection and background functions. Note that, for 
consistency with the Class code, the velocity source function 5e(r, fc) is given by 0(fc) = kV(k), where V(k) is the 
velocity perturbation in the Newtonian gauge. Primes indicate derivatives with respect to the argument. The index 
i refers to the redshift bin around reference redshift 2 ,> For the integrated terms “Len”, “G4” and “G5”, we have 
introduced 


Wt(r) 

W z G 4 (t) 

W G 5 (r) 


dfWiif) 


dfWi(f)k 


dfWi{f) 


2-5 s 

2 ) ( T o - t)(t 0 - f) 

2-5 s 

tq-t 

H' 2-5 s r 

1+ ^IP + (T 0 -r)aH +5S ~ U 


(A8) 


The window function W,(t) is given by the product of the number of galaxies per solid angle and redshift dN/dz/dfl 
multiplied, in the spectroscopic case, by a tophat centered around z t and of the width of the bin and, in the photometric 
case, by a Gaussian function with standard deviation equal to the half-width of the bin (see e.g. [58]). For numerical 
convenience we avoided to express SV in terms of the time derivative that requires numerical derivatives of Einstein 
equations, instead can be obtained analytically from Einstein equations [94] . The corresponding expressions 
presented in [30] are recovered integrating by parts A^ 5i (neglecting boundary terms since they vanish as r —> 0 and 
are unobservable for r = ro) and redefining consistently A^ 1; and A G2i so that the sum of these transfer functions 
coincides with the form given in previous work and in ClassGal 1 . Note that the “G4” term contains the Shapiro 
time delay, while “G5” is related to the integrated Sachs-Wolf effect (using integration by parts, the derivative acting 
on the Bessel function can be interchanged with the time derivative of the sum of Bardeen potentials). 


2. Modifications of CLASS 


We now describe how to implement the /nl parameter (and its running 71ng) for local primordial n on-G aussianity 
(PNG). In Class, it is convenient to use the Poisson equation $ = — (3f2om^o/ 2a ^ 2 ) ^co i n Equation (20), to obtain: 


Ab(z, k) = — 2(6 g (z) — 1)/nl ( w- 

. •£* 


9(0) 2U 

9(z dec) 


(A9) 


for more details about the derivation of this result see [55]. Here /nl, «ng, and 5 ec are constants describing the 
PNG parameter, its running, its pivot scale and the critical density for ellipsoidal collapse, respectively. The linear 
and scale-independent galaxy bias in absence of PNG is given by the function bc(z). The variable 7Z = —3$ p /2 (here 
‘Fp is the primordial Bardeen potential), is the primordial potential of curvature perturbations. The linear growth 
factor g(z) = (1 + z)D(z), where D(z) is the amplitude of the growing mode, appears since CLASS, by normalizing all 
source functions with respect to TZ (e.g., S(z,k)/TZ(k)), would otherwise evaluate $ of Equation (19) at decoupling, 
while in LSS convention it is evaluated at z = 0. In other words, we take into consideration the following relation 
between LSS and CMB conventions |5S] [HD] : 


f LSS 

/nl 


g(zdec) ,CMB 

5(0) /NL ’ 


(A10) 


and we use /nl = /nl S - 


http://cosmology.unige.ch/content/classgal 
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We refer to CLASS V2.3.2 2 , where a constant galaxy bias is implemented in the transfer. c module as a rescaling of 
the density transfer function. As a first step, we remove this constant from the transfer . c module. As described 
in [SO]: the code computes the density perturbation in the comoving gauge, related to the longitudinal Newtonian 
gauge by: 


Sen = 


long 


A 

J.,2 lon S > 


(All) 


where $i ong = kV\ ong is the divergence of the velocity in the longitudinal gauge. Furthermore, i5 co is a gauge-invariant 
variable and at first order in perturbation theory it coincides with the density perturbation in the synchronous-gauge 
comoving with dark matter used in |.‘I8j . Then, we modify the perturbations . c module to rescale the matter density 
source function S co (z, k) /7 Z(k) by a function &g(~) specified within the same module. Finally, we add the quantity 
A b ■ S co /n from Equation ( |A9| to the rescaled density contrast bc,(z)8 co /H. 

With this modifications, and allowing also for a time-dependent magnification bias s(z), the transfer functions given 
in Equation ([2]) become the terms given in Equation (A71. 
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